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Abstract. Aims. Stellar parameters and abundances of Na, Mg, Al, K, Ca, Sr, Ba, and Eu are determined for four very metal- 
Kf I poor (VMP) stars (-2.66 < [Fe/H] < -2.15). For two of them, HD 84937 and HD 122563, the fraction of the odd isotopes of 

Ba derived for the first time. 

Methods. Determination of an effective temperature, surface gravity, and element abundances was based on non-local thermo- 
dynamic equilibrium (non-LTE) line formation and analysis of high-resolution (R ~ 60000 and 90000) high signal-to-noise 
(S/N > 200) observed spectra. A model atom for H i is presented. An effective temperature was obtained from the Balmer H(, 
and H^ line wing fits. The surface gravity was calculated from the Hipparcos parallax if available and the non-LTE ionization 
balance between Cai and Can. Based on the hyperfine structure affecting the Ban resonance line /i4554, the fractional abun- 
dance of the odd isotopes of Ba was derived from a requirement that Ba abundances from the resonance line and subordinate 
lines of Ban must be equal. 

Results. For each star, non-LTE leads to a consistency of T^if from two Balmer lines and to a higher temperature compared to the 
(s^ ■ LTE case, by up to 60 K. Non-LTE effects are important in spectroscopic determination of surface gravity from the ionization 

I^N , balance between Ca i and Ca n. For each star with a known trigonometric surface gravity, non-LTE abundances from the lines of 

jrt ' two ionization stages, Ca i and Can, agree within the error bars, while a difference in the LTE abundances consists of 0.23 dex to 

0.40 dex for different stars. Departures from LTE are found to be significant for all investigated atoms, and they strongly depend 
on stellar parameters. For HD 84937, the Eu/Ba ratio is consistent with the relative solar system r-process abundances, and 
the fraction of the odd isotopes of Ba, fojj, equals 0.43±0.14. The latter can serve as an observational constraint on r-process 
models. The lower Eu/Ba ratio and fojj = 0.22±0.15 found for HD 122563 suggest that the s-process or the unknown process 
has contributed significantly to the Ba abundance in this star. 

Key words. Line: formation - Line: profiles - Stars: abundances - Stars: fundamental parameters - Stars: late-type - Stars: 
individual: HD 84937: individual: HD 122563: individual: BD-l-3°740: individual: BD-13°3442 

1 . Introduction early generation of stars. Determination of element abundances 

m in VMP stars is important for understanding nucleosynthesis 

Very metal-poor stars (hereafter VMP stars, with [Fe/HV] < -2) . ^, i r^ i t-. » i 

■' t- V > L / j_i y processes m the early Galaxy. The more metal-poor stars are, 

belonging to the halo population of the Galaxy were formed j- . . j .u jic i. . 

° ° ^ ^ •' on average, more distant, and they are more dimcult to ana- 

at the earliest epoch of star formation. These stars preserve , ^. u . -.u . i u j i . i 

^ ^ lyze than nearby stars with metal abundances close to solar, 

in the atmosphere the chemical composition produced by the , .• i ^u .u j r r •. j . • .• 

^ ^ ^ -^ In particular, the methods of surface gravity determination are 

T , rr ~ T~ r \ii u 1 • of limited usefulness. For distant field halo stars, trigonometric 
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wings, and, therefore, they are not sensitive to surface gravity. 
In fact, surface gravities of VMP stars can only be found from 
the ionization balance between neutral and ionized species of a 
selected atom. One mostly selects iron and Fe peak elements. 
In two recently found hyper metal-poor stars (hereafter HMP 
stars, with [Fe/H] < -5), only Ca is observed in two ionization 
stages (Christlieb et al. 2002, Frebel et al. 2005) and, therefore, 
can be a potential tool for the derivation of log g. 

The atmospheres of VMP stars are characterized by low 
electron number densities and low opacities, which are partic- 
ularly low in the ultraviolet spectral range. Hence, local ther- 
modynamical equilibrium (LTE) is not fulfilled in the atmo- 
spheric layers where the lines are formed. As was proven in 
the past decade, the surface gravity based on the ionization bal- 
ance Fei/Feii (Thevenin & Idiart 1999 Korn et al. 12003 ) and 
element abundances (for review, see Asplund [20051 1 are signif- 
icantly affected by departures from LTE. In this study, we se- 
lect four metal-poor stars observed with the High Dispersion 
Spectrograph of the Subaru Telescope and perform detailed 
analysis of stellar parameters and element abundances based on 
the non-local thermodynamic equilibrium (non-LTE) line for- 
mation for 10 important chemical species, including Hi, Cai, 
Can, etc. One of the stars, HD 122563, is a giant, and the re- 
maining three, HD 84937, BD -h3°740, and BD -13°3442, are 
at the hot end of the stars that evolve on time scales compara- 
ble to the Galaxy lifetime. Similar intrinsically bright objects 
are mainly found among VMP stars due to the selection effect. 
For example, one of two HMP stars, the Christlieb's star is a 
giant with \ogg = 2.2 (Christlieb et al. |2004| l and another one 
has an effective temperature T^ff = 6180 K (Aoki et al. 2006). 
In the selected stars, it is possible to study the formation of 
some of the strongest spectral lines, such as Can /18498, Mgi 
/15172 and /15183, which are among the few spectral lines ob- 
served in HMP stars. At [Fe/H] < -2, they are presumably of 
purely photospheric origin. Thus, the stars of our small sample 
can serve as reference stars in studies of VMP and HMP stars. 

The paper is organized as follows. Observations of the se- 
lected stars are described in Sect.|2] Non-LTE calculations are 
found in Sect. [3] The solar element abundances necessary for 
further differential abundance analysis are derived in Sect.|4] 
from the selected spectral lines. In Sect.|5] we determine ef- 
fective temperatures from the wings of Balmer lines and check 
surface gravities from the ionization balance Cai/Caii. Non- 
LTE abundances of 8 chemical elements are derived in Sect.|6] 
For two stars, HD 84937 and HD 122563, both the Ban reso- 
nance and subordinate lines are available, and we determine in 
Sect. l7.2l the fraction of the odd isotopes of barium that gives an 
independent estimate of relative contribution of the r-process 
and main i-process to heavy element abundances in the star. 
Our recommendations and conclusions are given in Sect.|8] 

2. Observations and data reductions 

The spectroscopic observations for our program stars were car- 
ried out on Feb. 2, 2002 with the High Dispersion Spectrograph 
(HDS; Noguchi et al. |2002| l at the Nasmyth focus of the Subaru 
8.2m telescope. Signal-to-noise (S/N) ratio measured per pixel 
(0.9km/s) is 200 to 400 for the brightest star HD 122563 and 



is lower, down to 100, for the remaining three stars. We use 
also high-quality observed spectra covering Can 3933, Can 
8498, the Hnes of Al i and K i from the ESO UVESPOP survey 
(Bagnulo et al. 2005) for HD 84937 and HD 122583 and from 
the ESO/UVES archive (Program ID: 67.D-0106 by P Nissen; 
67.D-0439 and 68.D-0094 by F Primas) for BD H-3°740 and 
BD-13°3442. The Balmer lines H„ and U/j in HD 84937 
and BD -i-3°740 are studied using observational data obtained 
by Klaus Fuhrmann with the fiber optics Cassegrain echelle 
spectrograph FOCES at the 2.2m telescope of the Calar Alto 
Observatory in 1999 and 1995. The high quahty of the FOCES 
data reduction accuracy has been emphasized by Korn (i2002l l 
who presents evidence for his claim that systematic errors of 
the effective temperatures derived from the Balmer line wing 
fits due to data reduction (rectification and flatfielding) are be- 
low 50 K. This holds for HD 84937. For the BD +3°740 spec- 
trum, the error may reach a slightly higher value (~70 K), be- 
cause before 1997 all spectra were degraded by mode transfer 
noise in the optical fibers (Grupp 2003 ) that did not allow an 
S/N > 150. ..200. Characteristics of observed spectra are sum- 
marized in Table[T] 

The Subaru and ESO/UVES spectra were reduced with 
a standard MIDAS/ECHELLE package for order identifica- 
tion, background subtraction, flat-fielding, order extraction, 
and wavelength calibration. Bias, dark current, and scatted 
light correction are included in the background subtraction. 
The spectrum was then normalized by a continuum function 
determined by fitting a spline curve to a set of pre-selected 
continuum windows estimated from the Kitt Peak Solar Atlas 
(Kuruczetal. il984D . 

3. Non-LTE calculations 

All the investigated elements are assumed to be trace elements, 
which means that we obtain statistical equilibrium (SE) popula- 
tions for each of them while keeping the atmospheric structure 
fixed. In the stellar parameter range we are concerned with, 
such approach is justified both for metals and for hydrogen. 
Hydrogen is not an important continuous opacity source com- 
pared to the H" ions at temperature below 7000 K. Although 
it becomes a dominant electron donor in the atmospheres of 
metal-poor stars ([Fe/H] < -1), our non-LTE calculations show 
that the ionization balance between H n and H i is close to ther- 
modynamic equilibrium (TE) in the layers where the Balmer 
line wings form. For instance, the population ratio H n/H i de- 
viates by no more than 1 % from the TE value below log T5000 - 
-2.5 in the model with T^g = 6390K, logg = 3.88, and [Fe/H] 
= -2.66. Based on the self-consistent non-LTE modelling of 
a solar-type star. Short & Hauschildt ( I2005I I conclude that the 
non-LTE effects of all the light metals (including Na, Mg, Al, 
K, and Ca among others) on the model structure and flux dis- 
tribution are small. 

Calculations are performed with the plane-parallel, homo- 
geneous, LTE, and blanketed model atmospheres computed 
for individual stellar parameters using the code MAFAGS 
(Fuhrmann et al. 11997b . For the program metal-poor stars, a- 
enhancement is accounted for with the abundances of Mg and 
Ca determined in this study, the Si abundance that follows the 
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Table 1. Characteristics of observed spectra. 



Telescope / 


Spectral range (A) 


R 


S/N 


Objects 




spectrograph 












Subaru / HDS 


4120 - 5430, 5520 - 6860 


90000 


>200 


HD 122563 










> 100 


HD 84937, BD+3°740,BD- 


-13° 3442 


VLT2/UVES 


3850 - 8550 


80000 


>200 


HD 122563, HD 84937 




VLT2/UVES 


3730 - 4990, 4760 - 6840, 
6600 - 10600 


48000, 
60000 


>200 


BD+3°740,BD -13-3442 




2.2-m/FOCES 


4500 - 6700 


60000 


> 100 


HD 84937 








35000 


> 100 


BD +3°740 





Mg one, and the oxygen abundance taken from the literature. 
We assume [O/Fe] - 0.5 where no other data is available. It is 
worth noting that oxygen in cool stellar atmospheres plays a 
minor role as a donator of free electrons and as opacity source, 
and the uncertainty of its abundance does not affect the calcu- 
lated atmospheric structure. 

In order to solve the coupled radiative transfer and statis- 
tical equilibrium equations for metals, we use a revised ver- 
sion of the DETAIL program (Butler & Giddings [1985]) based 
on the accelerated lambda iteration, which follows the efficient 
method described by Rybicki & Hummer ( II99IIII992] |. Non- 
LTE computations for H i are based on the complete hneariza- 
tion method as described by Auer & Heasley (1976). We apply 
the code NONLTE3 originally treated by SakhibuUin Cl983l) 
and advanced later (Kamp et al. 120031) . The departure coeffi- 
cients are then used to compute the synthetic line profiles. The 
metal line list has been extracted from Kurucz' d 1 9921 ) compi- 
lation. 
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Fig. 1. Departure coefficients logfe, for the first seven levels of 
H I in the model atmosphere of BDh-3°740. 



3.1. Non-LTE line formation for H 1 
3.1.1. Model atom of Hi 

We consider non-LTE line formation for H i in order to pro- 
vide diagnostics for the Balmer lines H„ and H^ that are used 
for an accurate temperature determination in cool stars. The 
model atom of Hi includes levels with principal quantum 
numbers up to n < 19 and energies adopted from Wiese et 
al. ( II9661 I. Transition probabilities are taken from the Vienna 
Atomic Line Data base (Kupka et al. 1 1 9991 ) and if not avail- 
able they are computed using the approximate formula of 
Bethe & Salpeter (1957) as implemented by Johnson {1922J- 
Photoionization cross-sections are evaluated using the exact ex- 
pression for the hydrogen atom. Collisional rates include con- 
tributions from inelastic collisions with electrons and hydrogen 
atoms. For electron-impact excitation, the R-raatiix calcula- 
tions of Przybilla & Butler d2004al ) are used for the transitions 
between the energy levels with n < 1 and the approximation 
formula of Johnson d 19721 ) for the remainder. For electron- 
impact ionization, we apply the Seaton formula as described 
by Mihalas (1978). The data on inelastic collisions with neu- 
tral H particles available in literature for H i have been recently 
investigated by Barklem (12007b . He concludes that they all are 
of questionable quality. Here we use the Drawin's (.1968 ) for- 
mula as implemented by Steenbock & Holweger ( |I9841 l. Since 



it provides only an order of magnitude estimate, we constrain 
the efficiency of hydrogenic collisions empirically. The rates 
calculated using this formula were multiplied by a scaling fac- 
tor 5 h = (no hydrogenic collisions), 0.1, 1, and 2 in order to 
make the effective temperatures derived from the two Balmer 
lines Wa and H^ in the selected stars consistent. Our best esti- 
mate 5 h = 2 agrees with that of Przybilla & Butler (.2004b J . 

Radiation transfer is treated explicitly for all Balmer and 
Paschen lines, for the first seven Lyman lines, for every bound- 
bound transition with niow =4-8 and n„p < 1 3, and for the four 
continua including Lyman, Balmer, Paschen and Brackett ones, 
in total, for 70 transitions. Radiative rates for the remaining 
transitions are evaluated using LTE mean intensities obtained 
from the formal solution of the radiation transfer equation. For 
the transitions arising between the energy levels with n < 8, 
the line absorption profile is computed in detail and includes 
Stark broadening, radiative damping, self-broadening, Doppler 
broadening (both thermal and turbulent), and fine structure as 
implemented by Barklem & Piskunov (.2003) in their subrou- 
tine HLINOP. For the remaining transitions, Doppler profiles 
are assumed. 



3.1 .2. Departures from LTE for H 1 

Our calculations show that the behavior of the departure coeffi- 
cients bi - n^^^^lnY^ and the line source functions are similar 
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for all the models representing the atmospheres of the Sun and 
selected stars. Here, n^^^^ and n^^^ are the statistical equilib- 
rium and thermal (Saha-Boltzmann) number densities, respec- 
tively. In Fig.[T] the departure coefficients of the lowest 7 levels 
of H I are shown as a function of continuum optical depth tsoqo 
refeiTing to A - 5000A in the model atmosphere of BD+3°740. 
The ground state and the first excited level (n - 2) keep their 
TE level populations throughout the atmosphere. Departures 
from LTE for the n - 3 level are controlled by Ha- In the 
layers, where the continuum optical depth drops below unity, 
Hq. serves as the pumping transition resulting in a slight over- 
population (up to ^ 5%) of the upper level. But starting from 
log T5000 - -3 and further out, photon escape from Ha wings 
causes an underpopulation of the n - 3 level. The levels with 
n > 4 stay in nearly detailed balance relative to each other 

For the Balmer lines of particular interest, H„ and H^, non- 
LTE leads to slight weakening of the core-to-wing transition 
compared to the LTE case and significant strengthening of the 
core. The core-to-wing transition is of particular importance 
for temperature determinations, because this part of profile is 
sensitive to Teff variations. 

We find that classical theoretical model atmospheres fail 
to reproduce a half-width of the H„ and Hg lines in the Sun 
(Fig.|2] for Ha) and cool stars (Fig. [3] and |5]). Similar conclu- 
sion is drawn by Przybilla & Butler ( 2004b ) and Fuhrmeister 
et al. (I2OO6I 1 for Ha in the Sun. Therefore the line core within 
0.7A to 1 .5A, with the exact number depending on stellar pa- 
rameters, is not included in the fit. In the entire range of stellar 
parameters we are concerned with, we obtain weak non-LTE 
effects for the Hp profile beyond the core independent of the 
applied scaling factor Sh- For example, for the model repre- 
senting the atmosphere of BD +3°740 (T^eff = 6340 K, log^ = 
3.90, [Fe/H] = -2.65), a difference in the LTE and non-LTE 
(Sh =0.1) profiles equals 0.1% of the continuum flux at AA = 
1.7A from the line center and rapidly decreases with increas- 
ing AA. In fact, non-LTE effects can be neglected deriving the 
effective temperature from H/j. For Ha, the coiTesponding dif- 
ference equals 2.7% of the continuum flux at AA - 1.7A and 
becomes smaller than 0.1% only beyond AA = 6. 7 A. A differ- 
ence in the non-LTE H„ profiles computed with S'h = 1 and 2 
is smaller than 0.1% of the continuum flux everywhere beyond 
AA - 1 .7 A. Both profiles are 0.7% stronger compared to that 
computed with ^h = 0.1 at A/l = 1.7 A. The three profiles agree 
within 0.1% beyond AA = 3.2A. In Sect. l5.3l we inspect an in- 
fluence of the Sh variation on the derived Teff of BD +3°740. 
Similar effects are expected for other program stars. 



3.1 .3. Hq, in the solar spectrum 

Hereafter, the theoretical profiles of Hq, and H^g are computed 
as accurately as possible, including a convolution of the pro- 
files resulting from the different broadening mechanisms. The 
Stark broadening profiles from Vidal et al. (I1970II1973T | are em- 
ployed. For self-broadening, we use two recipes. The first one 
is based on the Ali & Griem ( .1966 J theory (hereafter AG) and 
the second one on the self-broadening formalism of Barklem 
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Fig. 2. Theoretical non-LTE flux profile (continuous line) of 
the Balmer Hq, line computed for Teff = 5780 K using the AG 
self -broadening recipe compared to the observed spectrum of 
the Kurucz et al. (1984) solar flux atlas (bold dots). The dif- 
ferences between observed and calculated spectra (O - C) are 
presented in the upper part of figure. 



et al. (2000' hereafter BPO). We apply the revised subroutine 
HYDLINE incorporated in the SIU code. 

Non-LTE modelling of Hq. has been tested by comparing 
with solar flux observations taken from the Kitt Peak Solar 
Atlas (Kurucz et al. 1984) . We use the grid of the MAFAGS 
model atmospheres with various Teff but with the fixed values 
of log g - AAA and solar metal abundance. The line core within 
+ 1.1 A is not included in the fit. We obtain smaller non-LTE ef- 
fects for Hq in the Sun compared to that found in our program 
metal-poor stars. The best non-LTE fit is achieved for T^f^iAG) 
- 5780 K when the AG recipe is applied (Fig.|2]i and for a lower 
temperature T^ffiBPO) = 5720 K in the case of the BPO recipe. 
In each case, the LTE calculated wings for the same tempera- 
ture are slightly deeper (by up to 0.3%) than the observed pro- 
file, indicating a 30 K lower effective temperature value. The 
Hj3 wings are strongly blended in the solar spectrum. They are 
satisfactorily fitted using Tf^ffiAG) = 5780 K. 

3.2. Non-LTE calculations for metals 

For metals, our present investigation is based on the non-LTE 
methods treated in our earlier studies and documented in a 
number of papers, where atomic data and the problems of line 
formation have been considered in detail. Table |2] lists the in- 
vestigated atoms and cites the related papers. For Can, the ad- 
vanced effective collision strengths are used from the recent 
/^-matrix calculations of Melendez et al. ( 120071 ). 

We note that the most accurate atomic data are used in our 
non-LTE calculations. For minority species such as Mgi, Alt, 
and Ca i, where departures from LTE are mainly caused by su- 
perthermal radiation of non-local origin below the thresholds of 
low excitation levels, photoionization cross-sections are from 
the Opacity Project calculations (Seaton et al. |1994| l. Their ac- 
curacy is estimated at the level of 10%. Such small uncertainty 
translates to the abundance error of no more than 0.01 dex. For 
the collision-dominated atoms Nai and Ki and the majority 
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Table 2. Atomic models used in this study where NL is the 
number of levels in the included species' . 



Atom NL Reference 

Sodium Nai 58 Baumuller et al. ( |1998t . 

Gehren et al. ( l2004t . 
Shi et al. ( l2004t 
Zhao et al. ( [1998) , 
Gehren et al. ( 1200^ 
Baumuller et al. (JT996) 
Zhang et al. ll2006a;2006b|l 
Mashonkina et al. (i2007ai > 



Magnesium 


Mgi 


85 


Aluminium 


Ah 


59 


Potassium 


Ki 


67 


Calcium 


Cai 


63 




Can 


37 


Strontium 


Srn 


40 


Barium 


Ban 


35 


Europium 


Eun 


632 



Belyakova & Mashonkina (1997), 
Mashonkina & Gehren (2001) 
Mashonkina et al. ( |1999t 
Mashonkina ('2000), 
Mashonkina & Gehren J2000b , 
Mashonkina & Vinogradova ( 12007 1 



' Each model atom is closed with the ground state of the next ioniza- 
tion stage, 

2 combined levels based on 176 energy levels of Eun known from 
laboratory measurements. 



species such as Can, Srn, Ban, and Eun, the main source of 
uncertainties is collisional data. In SE calculations, we account 
for inelastic collisions both with electrons and neutral H par- 
ticles. Hydrogen collisions are computed using the formula of 
Steenbock & Holweger ( 119841 ) with a scaling factor 5 h found 
empirically from their different influence on the diff'erent lines 
of a given atom in solar and stellar spectra. We apply 5h = 
(no hydrogenic collisions) for Srn, Ban, and Eun as esti- 
mated by Mashonkina & Gehren ( 120001 I200TT ) and Mashonkina 
& Vinogradova (20071), 5h = 0.002 for All and 5h = 0.05 
for Nai according to Gehren et al. (,2004), Sh - 0.05 for Ki 
(Zhang et al. I2006al I2006bl l, and 5 h = 0.1 for Mg i and Ca i/ii 
as recommended by Mashonkina et al. (I2007al) . 

4. Line list and solar abundances for a differential 
analysis 

The spectral lines used in our study are listed in Table[3] 
Hyperfine structure (HFS) and/or isotope structure (IS) is taken 
into account when necessary with the data from Moore ( 119721 
IS for Mgi ^4703, ^5528, and ^5711), Nortershauser et al. 
( [19981 IS for Can i8498), Mc William et al. ( [T9951 Srii 
i42I5), Mc William ( .1998. Ban /14554), and Lawler et al. 
( [200n Eun ^4129). For Mg, Ca, Sr, and Eu, we use the frac- 
tional isotope abundances corresponding to the solar system 
matter (Anders & Grevesse .1989) . For barium, the even-to-odd 
isotope abundance ratio was allowed to vary, and its value is 
determined in Sect. l7.2l for the two stars. 

Based on non-LTE line formation, we determine the solar 
element abundance from each individual line. The exception is 
the Fen lines, which are analyzed assuming LTE. The solar 
flux profiles from the Kurucz et al. ( .1984J Atlas are fitted using 
the fixed oscillator strengths and van der Waals damping con- 



stants collected in Table[3] The profile of the strongest Ca line 
Can/13933 cannot be fitted although the best available atomic 
data are used. We assume that its core and inner wings are influ- 
enced by the chromospheric temperature rise and non-thermal 
depth-dependent chromospheric velocity field that is not part 
of the MAFAGS model. For this line, we adopt the mean Ca 
abundance derived from the Cai lines. 

We find that the mean solar abundances derived in this 
study agree within the error bars with the corresponding me- 
teoritic abundances of Asplund et al. (2005). This gives con- 
fidence in the reliability of our results obtained in analysis of 
stellar spectra. Although analysis is made line-by-line differen- 
tially with respect to the Sun, it is not strictly differential. A line 
which is measurable in the selected metal-poor stars, is strong 
in the solar spectrum, and the derived solar element abundance 
depends not only on its g/-value but also on van der Waals 
damping parameter. 

5. Stellar parameters 

Our results both for stellar parameters and element abun- 
dances are based on line profile analysis. In order to compare 
with observations, computed synthetic profiles are convolved 
with a profile that combines instrumental broadening with a 
Gaussian profile and broadening by macroturbulence with a 
radial-tangential profile. For a given star, the macroturbulence 
value Vtnac was allowed to vary by ±0.5 km s"' (lo"). Our strat- 
egy in this section is as follows. 

1. Stellar effective temperature is determined from the hy- 
drogen H„ and H/? line wing fitting based on non-LTE 
line formation. As was shown by Fuhrmann et al. ( |I993l l, 
the Balmer line profiles are strongly temperature sensitive, 
while the variation with gravity and the metal abundance 
is rather small. Hq, in particular is very insensitive to g and 
[M/H]. This is the basic justification for using the Balmer 
lines as a temperature indicator 

2. Next, the surface gravity is derived from the Hipparcos 
parallax with a mass determined from the tracks of 
VandenBerg et al. ( [2000b . The exception is BD-13°3442 
with no Hipparcos parallax available. 

3. The ionization balance between Cai and Can is checked 
for the fixed temperature and current estimate of log g. We 
denote the mean abundance from the subordinate lines of 
Cai as the Cai abundance and from Can /13933 and Can 
/18498 as the Call abundance. If (Cai - Call) exceeds the 
abundance error, the gravity is revised. We note that the 
two Ca n lines behave oppositely with varying gravity. The 
resonance line has the van der Waals broadened wings 
even at [Ca/H] < -2, and the derived Ca abundance de- 
creases/increases with increasing/decreasing gravity. The 
IR line is strengthened with decreasing \ogg due to de- 
creasing the H" continuous absorption and due to amplified 
departures from LTE for Ca n. This gives an additional op- 
portunity to constrain stellar gravity from a comparison of 
the element abundances from the two Can lines. 

4. Simultaneously, the iron abundance is obtained from the 
Fe n lines assuming LTE, and the abundances of a-process 
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elements Mg and Ca are determined from non-LTE analysis 
of the Mg I andCai lines. 
5. A micro turbulence velocity y„iic is derived from the 
strongest Fe ii, Ca i, and Mg i lines sensitive at given stellar 
parameters to a variation in Vmic- 

If any of the obtained parameters, either Teft, or log g, or 
[Fe/H], or c-enhancement, deviates from the one adopted orig- 
inally, the model atmosphere is recalculated and the steps 1 
through 5 are repeated. Thus, stellar parameters are determined 
iteratively. We describe below how it is made for the individ- 
ual stars and order them by the level of the accuracy of stellar 
parameters T^ff and \ogg available in the literature, from the 
higher accuracy to the lower one. The final stellar parameters 
and element abundances are presented in Table|4] The mean el- 
ement abundance is given where more than one line is used in 
analysis, and the error bars quoted in Table|4] is computed as 
the standard deviation cr - -\j'L{x - Xi)^/(n - 1). 

5.). HD 84937 

This star with a bright magnitude of 8.28 is among the best 
studied halo stars. Its effective temperature was determined 
many times (for review, see Korn et al. 12003') both spectroscop- 
ically from Balmer line wing fitting and using the infrared flux 
method (IRFM) with surprisingly consistent results. According 
to the LTE analysis of Korn et al. ( I2003I I. T^ff (Balmer) = 
6346 K, while Alonso et al. (1 19961 1 and Melendez & Ramirez 
(l2004l l give reff(IRFM) = 6330 K and 6345 K, correspond- 
ingly. The HippARCos parallax, n = 12.44 ± 1.06 mas is known 
with an accuracy better than 10%. 

The Balmer line profiles are computed applying the two 
self -broadening recipes, AG and BPO. In non-LTE, we obtain 
Tes(AG) = 6380 K from H„ and 6350 K from H^. Based on 
a S /N ratio of the observed spectrum and a sensitivity of the 
Balmer lines to a variation of Teff , we estimate the uncertainty 
of Jeff arising from the profile fitting as 30 K for Hq, and 60 K 
for H^. The core within ±1 A for Hq. and ±0.7 A for H^ is 
not included in the fit. A mean temperature, Teff = 6365 K, is 
adopted as a final value. In the case of the self-broadening from 
Barklem et al. (|2000] |. a mean temperature T^s{BPO) = 6300 K 
is derived from two Balmer lines. The maximum error of T^s 
then is, probably, not larger than 70 K. 

A trigonometric gravity was carefully calculated by Korn et 
al. (120031 1. log g = 4.00 ±0.12. For a temperature Teif = 6365 K 
obtained in this study, log g increases by less than 0.01. For the 
model with Teff/log g - 6365/4.00, we find a good agreement 
between the non-LTE Cal abundance based on 14 lines of Ca i 
and the non-LTE Call abundance from two Ca ii lines: (Cal - 
Call) = 0.01. In contrary, an agreement between diff'erent lines 
is destroyed at the LTE assumption (see Table|4]i. 

The error of the abundance difference (Cal - Call) is mainly 
determined by the error of the Call abundance, because the er- 
ror of the Cal abundance is very small for the program stars. 
It constitutes a-j yfn where n - 15 to 10 for different stars, 
and cr does not exceed 0.06 dex (Table|4]i. For Can /i3933 in 
different stars, the abundance error is mainly due to the uncer- 
tainties of the van der Waals damping constant and the pro- 



file fitting, while the uncertainty of microturbulence velocity 
dominates the abundance error for Can /t8498. The van der 
Waals damping constants predicted by the perturbation theory 
of Barklem & O'Mara ( 1998) have the uncertainty A log Ce - 
0.05 - 0.18 dex according to Barklem & Aspelund-Johansson 
([2005] ) and Barklem ('2006^. For Can /13933 in the program 
stars, this translates to the abundance error of up to 0.06 dex. 

For Can /13933 in HD 84937, only the line wings are 
fitted because the observed Doppler core is very broad and 
most probably affected by the chromospheric temperature rise 
and non-thermal and depth-dependent chromospheric velocity 
fields similar to solar ones. The abundance error is estimated 
as 0.06 dex. For Can /J8498, AVmc = 0.2 kms^' leads to 
cr =0.04 dex. Thus, the differences of LTE abundances (Call 
(/18498) - Cal) and (Call (A8498) - Call (^3933)) exceed their 
three to two errors. 

The iron abundance [Fe/H] = -2.15 and microturbulence 
velocity Vmic = 1.6 km s 'are determined from 10 Feii lines 
requiring that the derived abundances do not depend on line 
strength. The a-enhancement of the model is defined by the 
non-LTE Mg abundance, [Mg/Fe] = 0.39. 



5.2. HD 122563 

HD 122563 is one of the best observed halo stars. The 
SIMBAD Astronomical Database finds more than 300 refer- 
ences to studies of this star since 1980. However, stellar pa- 
rameters derived by different authors show a large spread of 
data. Most papers rely on an infrared flux method temperature 
reft(IRFM) = 4572 K, as determined by Alonso et al. a 999] ). 
Fulbright (2000) gets the lower value, T^^ - 4425 K, requir- 
ing the same iron abundance from the Fe i lines with high and 
lower excitation potential. A mean temperature Teif = 46 1 5 K 
is deduced by Barbuy et al. (120031) from using the b-y, V - R, 
V - K, and J - K calibrations of Alonso et al. (119991) . Surface 
gravity derived in different studies, either from the Hipparcos 
parallax or from the ionization balance between Fe i and Fe ii, 
ranges between log^ =1.5 (Barbuy et al. 120031 ) and \ogg — 0.6 
(Fulbright 2000). Significant discrepancies are seen for the mi- 
croturbulence velocity: Vmic varies between 2 kms"' (Barbuy 
et al. 120031 ) and Vmic = 2.9 km s ' (Simmerer et al. l2004] ). 

In spectroscopic determination of Teff, we use the observed 
spectrum taken from the ESO UVESPOP survey (Bagnulo et 
al. 120051 ). The wings of Hq. and H^ in HD 122563 are rather 
weak as compared to the turnoff stars. However, they are still 
sensitive to Teff. This is illustrated in Fig. [3] Non-LTE leads 
to a consistency of the temperatures derived from the two 
Balmer lines. For the two broadening recipes, we find T^f^iAG) 
= 4720 K and T^ftiBPO) = 4600 K. The uncertainty of T^s 
arising from the profile fitting is estimated as 60 K for Hq, and 
50 K for Hjj based onaS /N ratio of the observed spectrum and 
a sensitivity of the Balmer lines to a variation of Teff. The line 
core within ±1.5 A for Hq, and within ±1.2 A for H^ is not in- 
cluded in the fit. Figure|3]shows the best fits and the differences 
between observed and calculated spectra, (O - C). At the LTE 
assumption, Hq^ gives a 60 K lower temperature. 
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Fig. 3. Theoretical non-LTE flux profiles computed using the 
AG (continuous line for Tgff = 4720 K, dashed line for Teff = 
4600 K) and BPO (dotted line, Tetj = 4600 K) self-broadening 
recipes compared to the observed UVES profiles (bold dots) of 
Hq. (top panel) and Hp (bottom panel) in HD 122563. The (O - 
C) values are shown in the upper part of each panel. 



Fig. 4. The best non-LTE (5h - 0.1, continuous line) and 
LTE (dotted line) fits of the observed profiles (bold dots) of 
Can /18498 in HD 122563 (top panel) and BD-i-3°740 (bottom 
panel). The obtained non-LTE and LTE abundances of Ca are 
presented in TableH] Note, that it is impossible to achieve well 
fitting in LTE. 



As a final value we adopt Teff - 4600 K that is supported 
both spectroscopically and photometrically. The uncertainty of 
120 K is estimated as the difference between Teff from the two 
broadening recipes, because a treatment of the self-broadening 
of the hydrogen lines turns out to be the most significant source 
of Teff errors. The effect of Teff variation by 120 K on metal 
abundances is shown in Table|5] We note that a significant dis- 
crepancy of 0.33 dex appears between Mgi /14571 and the re- 
maining five Mgi lines when Tf^g increases by 120 K. This is 
due to a larger increase of departures from LTE for the line aris- 
ing from the ground state compared to that for the subordinate 
lines. 

A trigonometric gravity based on the Hipparcos parallax, 
j: = 3.76 + 0.72 mas, was carefully calculated by Barbuy et al. 
( l2003T l using Teft = 4600 K, log^ = 1.53 ± 0.2. Based on non- 
LTE hne formation, we check the ionization balance between 
Ca I and Ca ii for three gravitis, log g - 1 .50, 1 .30, and 1 .70, in 
order to find a right value. In each case, we determine the iron 
abundance from 9 lines, the Mg abundance from 6 lines, the 
Cal abundance from 15 lines, and the microturbulence velocity 
from the strongest lines of Fe ii (multiplet 42), Mg i (5 lines), 
and Cai (3 lines). The Can resonance line is most probably 
affected by the chromospheric temperature rise, and the Call 



abundance of HD 122563 is obtained only from Can /18498 
which is believed to be of photospheric origin. Strong depar- 
tures from LTE are found for Can /18498, such that, with the 
LTE assumption, its profile cannot be fitted by a variation in 
the Ca abundance. The LTE abundance is determined from line 
wing fits as shown in Fig.|4] We note that the theoretical LTE 
profile in Fig.|4](top panel) has been computed for an 0.07 dex 
larger Ca abundance, compared to that for the non-LTE pro- 
file, and this only helps to reproduce the wings but not the core 
of Can /18498. If to use the equivalent width of /18498, LTE 
leads to a 0.18 dex overestimated Ca abundance. In non-LTE, 
the core of /i8498 in HD 122563 is better fitted though not per- 
fectly. The line is saturated, and its half-width is only weakly 
sensitive to a variation of the Ca abundance and Vmic- We rely 
therefore on fitting the observed profile beyond AA =0.4A. The 
uncertainty of the profile fitting translates to the abundance er- 
ror of 0.06 dex. The maximal error of the van der Waals damp- 
ing constant A log C(, = 0.2 leads to A log s - 0.06 dex. 

The effect of gravity variation on the abundance of Ca from 
two ionization stages is shown in Table|5] We find that the Cal 
and Call abundances are consistent within the error bars for 
each of three gravities. The abundance difference is minimal 
for log ^ = 1 .50, (Cal - Call) = -(-0.04 dex and equals -h0.07 dex 
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and -0.08 dex for log g - 1.30 and 1 .7, correspondingly. Thus, 
a trigonometric gravity is supported spectroscopically, and we 
adopt log g = 1 .50 + 0.2 as a final value. 

5.3. BD +3°740 

For BD +3°740, we note a surprisingly large discrepancy be- 
tween the infrared flux method temperatures derived by Alonso 
et al. ( fT996] l. reif(IRFM) = 61 10 K, and Melendez & Ramirez 
(12004 1). TeffCIRFM) = 6440 K. An effective temperature de- 
termined in different studies from different spectroscopic and 
photometric methods (for review, see Ivans et al. 120031 1 varies 
between T^g = 6000 K (Ivans et al. l2003T l and T^g = 6355 K 
(Carrettaetal. l2002] l. 

Based on non-LTE line formation for Hi with hydrogen 
collisions taken into account applying 5 h = 2, we obtain from 
both Balmer lines Teff (AG) = 6340 K and T^g{BPO) = 6260 K. 
The uncertainty of T^g arising from the profile fitting is es- 
timated as 70 K for H„ and 80 K for H^. The core within 
±1.2 A for Hq, and within ±0.7 A for H^ is not included in 
the fit. Figure|5] shows the best fits and the sensitivity of the 
Balmer line profiles to a 100 K variation in T^g- For the LTE 
calculations, nearly the same quality fit of Hq, is achieved at a 
60 K lower temperature, however, provided that the core within 
± 1 .7 A is not included in the fit. In the inner wings (A/l - 1 .2A 

- 1.7A), the LTE profile is too broad for T^g obtained from 
the far wings (see Fig.|5]l. When non-LTE calculations are per- 
formed with reduced hydrogenic collisions (S'e = 0. 1), a 100 K 
larger effective temperature is determined from the Hq. wings, 
beyond 2A from the line center, while the inner wings are ob- 
tained to be too narrow. Thus, the uncertainty of reff(HQ) due 
to various treatment of hydrogen line formation is estimated as 
(+100 K, -60 K). Having in mind that non-LTE effects are neg- 
ligible for Wj}, a consistency of the temperatures derived from 
two Balmer lines favors the non-LTE approach with hydrogen 
colhsions taken into account with Sn of 1 to 2. 

As a final value, T^g - 6340 K, is adopted with the un- 
certainty of lOOK. The effect of T^g variation on metal abun- 
dances is shown in Table|6] 

According to Fuhrmann ( 119981 1, a trigonometric gravity 
corresponding to the Hipparcos parallax, n - 7.80 ± 2.09 mas, 
is inferred as \ogg - 4.18, and it can lie between \ogg - 3.92 
and \ogg - A-.lil due to the relatively large error of the paral- 
lax. Asplund et al. (2006) use the mean value of the absolute 
visual magnitude derived from the Hipparcos parallax and from 
Stromgren photometry and obtain log g - 4.04. In both studies, 
Teff = 6260 K is used. We have checked the ionization balance 
between Cai and Can for \ogg - 4.06 and 3.90, applying T^g 

- 6340 K. Table|6]presents the effect of gravity varying on the 
obtained abundance of Ca. The Cal abundance is determined 
from 10 lines with the measurement error cr - 0.05 dex. For 
Can /18498, AVmic = 0.2 kms"' leads to the abundance er- 
ror of 0.06 dex, while the uncertainty of the profile fitting to 
0.03 dex. For the Ca ii resonance line, the total abundance error 
is estimated as 0.06 dex. 

For log^ - 3.90, we find a consistency within 0.04 dex 
of the Cal and Call abundances. For the higher gravity, an 
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Fig. 5. The best non-LTE fits of Hq, (top panel) and Wp (bot- 
tom panel) in the observed FOCES spectrum (bold dots) of 
BD+3°740 achieved using the AG (continuous line, T^g - 
6340 K) and BPO (dotted line, T^g = 6260 K) self -broadening 
recipes. For comparison, we show (O - C) also for the best LTE 
AG fit of Hq computed for T^g = 6280 K. Note that the LTE 
inner wings (A/i - 1.2A- 1.6A) are broader compared to the 
non-LTE ones although the lower effective temperature is used 
in the LTE calculations. In order to illustrate the sensitivity of 
the line profiles to T^g, the non-LTE AG profiles are shown also 
for Teff - 6240 K (dashed curves). 



agreement between the different lines of Ca is destroyed: Call 
(^8498) - Cal = 0.09 dex and Call (^8498 ) - Call (^3933) = 
0.11 dex. As a final value, we adopt \ogg - 3.90±0.15 which 
is the lower limit of a trigonometric gravity and which is best 
supported spectroscopically. 

5.4. BD- 13° 3442 

This is the faintest star of our small sample with no Hipparcos 
parallax measured. In the observed UVES spectrum of this 
star, only Ho- can be used to derive a temperature. The ob- 
tained spectroscopic temperature based either on the Ali & 
Griem ( |T966l l theory, T^g(AG) = 6390 K, or on the BPO recipe, 
Tf.g{BPO) - 6310 K, is lower compared to an infrared flux 
method temperature, reff(IRFM) = 6484 K, as determined by 
Melendez & Ramirez ( E004] l. As a final value, T^g = 6390 K 
is adopted with the uncertainty similar to that for BD +3°740, 
100 K. 
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Using the absolute visual magnitude from Stromgren 
photometry, Asplund et al. ( 120061 1 obtain logg = 3.86 for 
BD -13°3442. A 0.02 dex larger surface gravity is obtained for 
a slightly higher temperature, Tf^g - 6390 K, determined in this 
study. Assuming that the gravity error is similar to that found 
for BD -i-3°740, we check the ionization balance between Cai 
and Can for logg = 3.88 and 3.60. For each of two gravities, 
the Cal and Call abundances turn out to be consistent within 
the error bars. The effect of a gravity variation on the obtained 
abundances of Ca is shown in Table|6] We adopt logg - 3.88 as 
a final value. It is worth noting that the calculated atmospheric 
parameters of BD -13°3442 and BD +3°740 are very similar. 

It is interesting to inspect the position of the selected stars 
on evolutionary tracks using the obtained stellar parameters. 
The grid of stellar models of VandenBerg et al. (2000) is used. 
According to interpolation along lines of constant radii in a 
cube of masses, metal abundances and [a/Fe] ratios for avail- 
able tracks we derive the HRD position of HD 84937 with a 
mass of 0.784 Mq and an age of 14.0 Gyr (linear interpola- 
tion) or an age of 13.7 Gyr (logarithmic interpolation). The 
age determined for the known trigonometric parallax is slightly 
higher, i.e. 15.3 Gyr. The other two turnoff" stars, BD+3°740 
and BD-13°3442, are outside the set of evolutionary tracks 
with their low metal abundances. Their masses and ages can 
only be estimated by extrapolation (the lowest metal abundance 
of the grid is at [Fe/H] = -2.3). Their masses are 0.766 Mq 
(BD-h3°740) and 0.790 Mq (BD-13°3442), whereas their ages 
are 14.9 Gyr (BD-h3°740) and 13.0 Gyr (BD-13°3442). The 
estimated ages are acceptable in view of the WMAP results 
that center around an age of the universe of 13.7 Gyr (Spergel 
et al. 120031) . HD 122563 lies on the giant branch. In such case, 
extrapolation leads to too uncertain results for stellar ages. 
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Fig. 6. Trends of non-LTE (filled circles) and LTE (open cir- 
cles) abundance with line strength determined from the Cai 
fines in HD 122563. Note the steep trend of the LTE [Ca/Fe] 
values with line strength above Wa - 60 mA. The mean non- 
LTE value is shown by dotted line. 



We find that HD 122563 reveals significantly lower ele- 
ment/Fe abundance ratios compared to the remaining stars. For 
the most investigated elements, there is underabundance by 
1 .2 dex for Eu, by more than 0.75 dex for Ba, 0.5 dex for Al, 
0.3 dex for K and Sr, and 0.2 dex for Ca. This reflects most 
probably a chemical inhomogeneity of the interstellar medium 
in the early Galaxy. An environment, where HD 122563 was 
formed, was enriched with the neutron-capture elements and 
the elements with odd nuclear charge Al and K to the less ex- 
tent compared to the interstellar gas out of which the remaining 
stars were formed. 



6. Abundance analysis 

In order to derive aluminium and potassium abundances, we 
use the UVES archive spectra because the All /13961 and Ki 
AA1665, 7699 lines are not covered by the available Subaru 
spectra. For BD +3°740, we cannot calculate Al abundance 
from the line at /13961A because the observed spectrum in 
this region is missing. The Ki A7699 line in BD-13°3442 is 
blended with telluric lines and therefore is excluded from anal- 
ysis. In BD -h3°740 and BD-13°3442, the Euii ^4129 and Ban 
/16496 lines cannot be extracted from noise. The Ba abundance 
of these stars is obtained from a single line, Ba ii /14554, as- 
suming a pure r-process even-to-odd Ba isotope ratio of 54 : 
46 (Arlandini et al. |1999l l. It is worth noting that the Ban line 
in these stars is weak (Wa ~ 10 mA) and is insensitive to a 
variation of Ba isotopic abundances. The derived non-LTE and 
LTE abundances are presented in Table|4] 

For the selected elements, non-LTE abundances of the 
two stars were determined in our earlier studies based on the 
FOCES observed spectra (Gehren et al. 120061 Na and Mg in 
HD 84937 and BD-h3°740; Mashonkina et al. 12007a. Mg and 
Ca in HD 84937). Shimanskaya & Mashonkina (l200TT l cal- 
culated the non-LTE abundances of Mg for HD 84937 and 
BD-i-3°740 based on the equivalent widths from the literature. 
In this paper, we obtain very similar results. 



6.1. Departures from LTE 

Significant non-LTE efifects are found for almost all lines in 
our program stars. Non-LTE abundance corrections Anlte - 
log enon-LTE -log fiLTE are negative for alkafi atoms Na i and K i. 
They grow in absolute value with the strength of the observed 
line, i.e. with decreasing effective temperature. The maximal 
values are reached for HD 122653: ANLTE(Nai) = -0.44 dex 
and Anlte(Ki) - -0.35 dex. 

For the photoionization dominated minority species Mgi, 
All, and Cai, Anlte is positive. We note a strong tempera- 
ture dependence of departures from LTE for All. Non-LTE 
abundance correction is equal 0.86 dex for the hottest star 
BD-13°3442, while consists of only 0.04 dex for the coolest 
star HD 122653. For the Mg i and Ca i lines, Anlte does not ex- 
ceed 0.2 dex. However, an advantage of the non-LTE approach 
is clearly seen in Fig.|6]where we show the Ca abundances cal- 
culated from the neutral Ca lines in HD 122653. Non-LTE re- 
moves the steep trend of the [Ca/Fe] values with line strength 
obtained under the LTE assumption above Wa = 60 mA. 

As was found theoretically in our earlier studies, non-LTE 
may lead either to strengthening or to weakening the Ban 
(Mashonkina et al. 1999) and Srii (Mashonkina & Gehren 
1200 11 1 lines depending on stellar parameters and element abun- 
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dance. This is nicely illustrated with our small sample of stars. 
A slightly negative non-LTE abundance correction is obtained 
for Srii /14215 in HD 122653 and a positive one for two the 
most metal -poor stars, while Anlte - for HD 84937. For the 
Ban lines in our program stars, Anlte is positive and reaches 
0.50 dex for BD-13°3442. Non-LTE effects lead to Anlte = 
-1-0.16 dex for Eu II /14129. 

6.2. Comparisons with non-LTE studies 

For the elements under unvestigation, systematic analyses 
based on non-LTE line formation are few in number. We find 
that the non-LTE abundances obtained in this paper agree 
within the error bars with the corresponding abundances in 
the literature for common stars. For the three stars, HD 84937, 
BD-h3°740, and BD-13°3442, the absolute abundance differ- 
ence (this study - Idiart & Thevenin, I2000I I does not exceed 
0.08 dex for Mg and 0. 12 dex for Ca. The non-LTE abundances 
of Mg and Ca given by Gratton et al. C2003J for HD 84937 are 
consistent within 0.06 dex with ours. 

The non-LTE abundance of Na determined by Mishenina et 
al. (I2003J for HD 84937 agrees within the error bars with our 
data. For HD 122563, we derive a 0. 14 dex larger Na abundance 
compared to that given by Andrievsky et al. (|20071 l. This is 
most probably due to the lower surface gravity {logg = 1.1) 
used in their study and, therefore, the stronger non-LTE effects 
for the Na I lines. 

The non-LTE abundance of potassium obtained by Takeda 



et al. ( 120021 ) agrees within 0.03 dex with our value for 
HD 122563. 

The non-LTE abundances of Al and heavy elements Sr, Ba, 
and Eu are determined for the first time for the program stars. 

6.3. Galactic abundance trends: non-LTE vs. LTE 

In the view of significant departures from LTE found for the 
program stars, it is important to investigate what might change 
from the abundance trends determined from the larger samples 
of stars analyzed in LTE in the literature. 

The a-process elements, Mg and Ca. There are well es- 
tablished observational evidences for that, in old metal-poor 
([Fe/H] < -1) stars of the Galaxy, Mg and Ca are overabun- 
dant relative to iron (see e.g. recent determinations of Fulbright 
120001 Cayrel et al. l2004l Cohen et al. ,2004. Barklem et al. 
120051 they all computed under LTE). The different samples of 
stars reveal close together average abundance ratios for Mg/Fe 
and Ca/Fe and the small scatter of data at the level of the 
measurement errors. E.g. Cohen et al. ( I2004I I obtain the mean 
abundances [Mg/Fe] = 0.39+0.13 and [Ca/Fe] = 0.31+0.12 for 
the sample of the turnoff stars with -3.5 < [Fe/H] < -2, and 
Cayrel et al. (|2004] | find [Mg/Fe] = 0.27+0.13 and [Ca/Fe] = 
0.33+0.1 1 for the sample of the cool giants with -4 < [Fe/H] 
< -2.4. In the range of stellar parameters the studies are con- 
cerned with, non-LTE leads to larger abundances of Mg and Ca 
derived from the lines of neutral atoms. According to Zhao & 
Gehren (.2000 . Mg i), Idiart & Thevenin (2000 Mg i and Ca i), 
and Mashonkina et al. (I2007al Ca i), non-LTE abundance cor- 



rections may reach 0.2 dex depending on stellar parameters. 
For the Cohen et al. sample, they are expected to be simi- 
lar to that calculated for our program turnoff stars BD+3°740 
and BD-13°3442, Anlte = +0.04 dex for Mg and Anlte = 
+0.14 dex for Ca. Based on our calculations for HD 122563, 
we expect that the LTE abundances of Mg and Ca in the cool 
giants of the Johnson ( I2002I I and Cayrel et al. ( I2004I I samples 
may be larger in non-LTE by approximately 0.2 dex. We note 
that accounting for departures from LTE leads to nearly identi- 
cal overabundances of Mg and Ca relative to Fe in the turnoff 
stars of Cohen et al. and the cool giants of Cayrel et al. with 
the Ca/Mg abundance ratio close to solar value. This makes the 
conclusions of the cited papers about a similarly constant ratio 
between the yields of iron and the c-process elements Mg and 
Ca in the early Galaxy more solid. 

Aluminium. As for neutral magnesium and calcium, non- 
LTE leads to larger abundances of Al derived from the res- 
onance lines of All compared to those from LTE analysis. 
However, departures from LTE for Al i depend strongly on ef- 
fective temperature and metallicity. As was shown above, the 
metal-deficient turnoff stars are affected most strongly by Al 
abundance coiTections reaching +0.86 dex for BD-13°3442. 
When non-LTE effects are neglected. Galactic abundance 
trends have a large scatter of data for the stars of close metallic- 
ity. E.g. for their sample of VMP stars (-3.8 < [Fe/H] < -1.5), 
Barklem et al. ( 2005) find the scatter of the order of 0.5 dex for 
the [Al/Fe] and [Al/Mg] ratios, while the scatters are less than 
the typical abundance errors for [Mg/Fe] and [Ca/Fe]. The data 
of Fulbright (120021) on [Al/Fe] reveals even larger scatter up to 
1 dex at [Fe/H] < -2. 

Some abundance analyses (e.g. Carretta et al. 120021 Cayrel 
et al. 2004, Cohen et al. 2004) carried out under the LTE as- 
sumption use the non-LTE abundance corrections published 
in Baumiiller & Gehren (I1997I I. We note that Baumiiller & 
Gehren d 19971 ) give Anlte for the small grid of model atmo- 
spheres, and individual corrections obtained by interpolating 
for given stellar parameters should be considered as a coarse 
treatment of the non-LTE effects. Nevertheless even a rough 
approach leads to a dramatical change of Galactic abundance 
trends. For the sample of stars in the metallicity range -3.6 < 
[Fe/H] < -2.4, CaiTetta et al. obtain the average LTE abun- 
dance [Al/Fe] = -0.59+0.15, while the Al abundance increases 
up to [Al/Fe] = +0.10+0.22 when the non-LTE corrections are 
applied. A different case is the Cayrel et al. study. They use an 
inappropriate coiTection for their sample of cool giants, Anlte 
- +0.65 dex that corresponds to the model with Teif - 5500 K, 
log g = 3.5, [Fe/H] = -3. Our calculations for HD 122563 show 
that the Al abundance correction is much smaller, Anlte - 
+0.04 dex. The non-LTE abundances of Al presented by Cayrel 
et al. in their Fig. 8 are therefore overestimated. 

Alkalis, Nai and Ki. Contrary to Mg, Ca, and Al, non-LTE 
leads to smaller abundances of sodium and potassium com- 
pared to those determined in LTE. As was shown by Baumiiller 
et al. ( 1998) and Andrievsky et al. (,200 7 ) and supported in this 
paper, departures from LTE for Na i depend strongly on the sur- 
face gravity and the element abundance of a star. For the cool 
giant HD 122563, the non-LTE abundance correction of Nai 
/15889 / /15895 is a 0.15 dex / 0.20 dex more negative com- 
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pared to that for the turnofF star BD- 13° 3442 of close metal- 
licity and Na abundance. For our program turnofF stars, Anlte 
consists of -0.47 dex / -0.39 dex for the atmospheric parame- 
ters that match HD 84937 and decreases in absolute value with 
decreasing metallicity, by approximately 0.2 dex for A[Fe/H] 
= -0.5. As for Nai, the non-LTE effects for Ki increase with 
decreasing surface gravity. For alkalis, a full treatment of non- 
LTE line formation is therefore required for each individual 
star Andrievsky et al. (2007^ show a clear advantage of such 
an approach compared to the use of a constant non-LTE abun- 
dance correction. They achieve the small star-to-star variation 
in the [Na/Fe] and [Na/Mg] ratios comparable with the mea- 
surement errors and consistent results for the samples of VMP 
dwarfs and unmixed giants. 

Stellar abundance ratios between the odd-Z and even-Z ele- 
ments such as Na/Mg, Al/Mg, and K/Ca are important for test- 
ing current nucleosynthesis models. We emphasize that exactly 
these abundance ratios are most affected by departures from 
LTE. E.g. non-LTE leads to an increase of the Al/Mg ratio in 
the program turnoff stars by 0.77 dex. In contrary, the Na/Mg 
ratio decreases, on average, by 0.42 dex. It is worth noting that 
our non-LTE abundance ratios Al/Mg and Na/Mg agree well 
with the Galaxy chemical evolution calculations of Timmes et 
al. (119951 for the corresponding metallicity range. An excep- 
tion is a too low value of [Al/Mg] in HD 122563. As was noted, 
this star reveals an overall deficiency of the elements requiring 
for their synthesis a large flux of neutrons. 

Our data on potassium with the mean non-LTE ratio [K/Ca] 
= -0. 19+0.06 favor the Galaxy chemical evolution calculations 
of Goswami & Prantzos (120001) . Similar conclusion was drawn 
by Shimansky et al. (l2003l l and Zhang et al. (I2006bl l based 
on non-LTE determinations of the potassium abundance in the 
mildly metal-poor ([Fe/H] > -1) stars. 

Due to the strong influence of non-LTE effects in the at- 
mospheres of metal-deficient stars LTE abundance analyses of 
Na, Al, and K are nearly useless for an investigation of Galactic 
chemical evolution. 

Neutron-capture elements Sr, Ba, and Eu. The LTE abun- 
dance analyses (see e.g. Fulbright 120021 Carretta et al. ,2002, 
Barklem et al. 120051 ) determine a steep decline of the [Ba/Fe] 
abundance ratios at [Fe/H] < -2. Based on our non-LTE calcu- 
lations we expect that the trend might be shallower in non-LTE. 
E.g. for the program turnofF stars, [Ba/Fe] decreases from 0.00 
at [Fe/H] = -2.15 down to -0.23 at [Fe/H] = -2.66 in non-LTE, 
while from -0.15 down to -0.73 in LTE. 

As found in many studies (for a review of the literature 
under 2004, see Travaglio et al. 120041 and the later study of 
Barklem et al. 120051 they all computed under LTE), the abun- 
dance trends of Sr/Fe, Ba/Fe, Eu/Fe, and the ratios among 
neutron-capture elements have a large scatter at fixed metal- 
licity that appears to increase with decreasing [Fe/H]. A part 
of the dispersions may be due to ignoring deviations from LTE 
in abundance determinations, because the non-LTE effects for 
the lines of Srii and Ba ii depend strongly on stellar parameters 
and the element abundance. 

The measured scatter is usually interpreted as cosmic in ori- 
gin, because iron and neutron-capture elements are produced 
in different sites. We suggest to plot the ratios among neutron- 
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r-process abundance (Sr/Ba)^ - 0.62 according to Arlandini 
et al. ( 119991 1 by sohd fine. 



capture elements against the abundance of Ba, but not the Fe 
abundance. In Fig. |2l we combine the Sr/Ba abundance ra- 
tios from Barklem et al. (120051) with our data. Except for the 
few outliers, the scatter at fixed Ba abundance does not exceed 
0.5 dex at [Ba/H] > -3.5 and approaches to 1 dex at the lower 
[Ba/H]. This makes possible to see a growth of Sr/Ba towards 
the lower Ba abundance. For comparison, the scatter in the run 
[Sr/Ba] - [Fe/H] (see Fig. 22 from Barklem et al.) approaches to 
1 .5 dex at [Fe/H] < -2.5 and hides the correlations between the 
abundance ratios. According to the recent study of Mashonkina 
et al. (2007b), the Zr/Ba and Y/Ba abundance ratios in halo 
stars grow towards the lower Ba abundance. A clear abundance 
trend seen in Figure|2]indicates an existence of the process that 
yielded in the early Galaxy strontium with little contribution to 
Ba. The present study supports the ideas of Truran et al. (120021 ) 
and observational findings of Aoki et al. (120051) and Honda et 
al. (120061) that the light and heavy neutron-capture elements in 
VMP stars come from distinct processes. Synthesis of the light 
neutron-capture elements remains a challenge to theorists. 

7. The r/s-process controversy for HD 84937 and 
HD 122653 

7.1. [Eu/Ba] abundance ratios 

It is well known (among many observational studies see e.g. 
Spite & Spite [19781 and Barklem et al. l2005l ) that in old 
metal-poor galactic stars europium is overabundant relative to 
Ba. This is interpreted as a a dominance of the r-process in 
production of heavy elements beyond the iron group in the 
early Galaxy. The solar system abundance ratio of Eu to Ba 
contributed by the r-process relative to the total abundances, 
[Eu/Ba],-, equals 0.70 according to Arlandini et al. ( I19991 ). 
In several studies, Sneden et al. ( 119961 120001 ), Cowan et al. 
( fT999l 12002] ), and HiU et al. ( |2002] ) have presented arguments 
supporting constant relative r-process element abundances dur- 
ing the history of the Galaxy (at least, where 56 < Z < 70). 
In two stars, BD-h3°740 and BD-13°3442, flie Eu abundance 
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is not determined. Our computations for [Eu/Fe] - 0.5 (which 
translates to [Eu/Ba] = 0.7) predict the central depth of Euii 
A4129 in these stars at a level of 0.5%. Such a weak line cannot 
be measured in the spectra available to us. Thus, the Eu abun- 
dance in BD-h3°740 and BD-13°3442 may be large enough to 
give a [Eu/Ba] ratio close to 0.7. We obtain that HD 84937 re- 
veals abundances close to a pure r-process production of heavy 
elements with [Eu/Ba] - 0.70. For HD 122563, the correspond- 
ing abundance ratio is smaller, [Eu/Ba] = 0.50, that suggests a 
contribution from the main i-process to observed barium. 

7.2. The fraction of the odd isotopes of Ba 

A determination of the Ba even-to-odd isotope abundance ra- 
tio in stars becomes possible due to the significant hyperfine 
structure (HFS) affecting the Ban resonance lines of the odd 
isotopes. In total, the resonance line /14554 has 15 components 
spread over 58 mA. The larger the fraction of the odd isotopes, 
the stronger the HFS broadening of /14554 and the larger the 
energy absorbed in this line. The HFS broadening of the Ban 
subordinate lines is small. The maximal effect on the derived 
abundance is found for A6496 and it does not exceed 0.01 dex 
for the Sun and is much smaller for our program metal-poor 
stars. We use the method described in detail in our previous 
study (Mashonkina & Zhao , 2006 ) and find the Ba even-to-odd 
isotope abundance ratio in a given star from the requirement 
that Ba abundances derived from the Ba ii resonance line and 
the subordinate lines must be equal. 

For both stars, HD 84937 and HD 122563, the only subor- 
dinate line. Ban /i6496, is used to determine total Ba abun- 
dance. It is given in Table|4] The Ban /15853 line cannot be 
extracted from noise in HD 84937 and is affected by a hot pixel 
in HD 122563. We do not use Ban /16141 due to blending with 
a Fe I line. 

Having fixed the Ba abundance of HD 84937 ([Ba/Fe] = 0) 
and varying the even-to-odd isotope abundance ratio, we ob- 
tain a consistent element abundance from A4554 for a frac- 
tion of the odd isotopes fodd = 0.43 ± 0.14. The uncertainty 
of the desired value is mainly caused by the uncertainty of the 
microturbulence velocity, which is estimated as -H0.2kms"^ 
Recent r-process calculations of Kratz et al. (120071) predict 
fodd = 0.44. A close value, fodd = 0.46, is obtained from 
r-residuals in the stellar model of Arlandini et al. (1999). Thus, 
both the [Eu/Ba] abundance ratio and the even-to-odd Ba iso- 
tope abundance ratio suggest a pure r-process production of 
heavy neutron-capture elements in HD 84937. 

For HD 122563, we find f„dd = 0.22 ± 0.15 with the un- 
certainty mainly caused by the uncertainties of T^ff (120 K) 
and Vmic (0.1 kms"') (see Table|5]l. This fraction is close to the 
solar system value (0.18). Taking into account the Icr error, 
we estimate the r-process contribution to Ba in HD 122563 to 
be no larger than 70% in agreement with the [Eu/Ba] abun- 
dance ratio in this star. The only halo star with a measured 
fraction of the odd isotopes of Ba is HD 140283 (Lambert & 
Allende Prieto l2002b . Similarly to HD 122563, it shows a lower 
value fodd = 0.31 compared with that of two other halo stars 
with a measured fraction of the odd isotopes of Ba, HD 84937 



and HD 103095 (fodd = 0.42, Mashonkina & Zhao |2006] |. 
We note that both HD 140283 and HD 122563 reveal lower 
neutron-capture element abundances compared to those of the 
stars of close metallicity and look as r-process-poor stars. For 
HD 140283, [Ba/Fe] = -0.80 according to Mashonkina et al. 
(I2007b[ l. Using the same stellar parameters Teff - 5810 K, 
logg = 3.68, and [Fe/H] = -2.43, we estimate [Eu/Fe] = -0.31 
based on the UVESPOP observed spectrum (Bagnulo et al. 
120051) and non-LTE line formation. The obtained ratio [Eu/Ba] 
= 0.49 is smaller than the value typical of the halo stars (e.g. 
[Eu/Ba] - 0.69 deduced by McWiniam |1998l as the mean ratio 
for the sample of stars with [Fe/H] < -2.4) and smaller than 
the relative solar system r-process abundance [Eu/Ba]^ = 0.70 
(Arlandini etal. [T999] l. 



8. Conclusions 

In this paper, we determine for the first time stellar parameters 
Teff and logg of four VMP stars (-2.15 < [Fe/H] < -2.66) 
based on non-LTE line formation. Effective temperature is de- 
rived from the Balmer Hq, and H^g line wing fits. For each star, 
non-LTE leads to a consistency of T^ff from two lines and to a 
higher temperature compared to the LTE case, by up to 60 K. 

We find a clear advantage of the non-LTE approach in spec- 
troscopic determination of surface gravity compared to the LTE 
case. For each of three stars with Hipparcos parallax available, 
a surface gravity obtained from the non-LTE ionization bal- 
ance between Ca i and Ca ii agrees within the error bars with 
a trigonometric one. The Cai/Caii line intensity ratio involv- 
ing Ca II /18498 is particularly sensitive to a variation of sur- 
face gravity. However, the wings of the Ca ii resonance lines 
can also be used for spectroscopic determination of log g of 
VMP stars. We conclude that non-LTE calculations based on 
ID model atmospheres provide reliable results for Teff and 
logg of VMP stars in the wide range of temperatures (Teff - 
4600 K to 6400 K) and gravities (logg = 4.0 down to 1.5). 

Due to strong influence of non-LTE effects in the atmo- 
spheres of VMP stars abundances of Na, Mg, Al, K, Ca, Sr, 
Ba, and Eu cannot be determined in a reliable way assuming 
LTE. Departures from LTE depend strongly on stellar parame- 
ters such that a precise abundance analysis is only possible on 
the base of non-LTE calculations for each individual star. 

We derive the fraction of the odd isotopes of B a in two stars, 
HD 84937 and HD 122563, and, thus, increase the number of 
halo stars with a measured fodd to four. Two stars among four 
have low fodd- HD 122563 with fodd = 0.22+0.15 (this study) 
and HD 140283 with fodd = 0.31 (Lambert & Allende Prieto 
120021 ). Both stars reveal the lower Eu/Ba ratios compared to 
the relative solar system r-process abundances. These results 
suggest that the i-process has significantly contributed to the 
Ba abundance in these stars, or that the unknown process that 
produced high Sr/Ba ratios in these stars yielded the low fodd 
value. Two remaining stars have the larger fraction of the odd 
isotopes of Ba: fodd = 0.43+0.14 in HD 84937 (this study) and 
fodd = 0.42 in HD 103095 (Mashonkina & Zhao|2006|. In both 
stars, the Eu/Ba ratios are consistent with the relative solar sys- 
tem r-process abundances. We conclude therefore, that fodd 
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found in these stars can serve as observational constraint on 
r-process models. 
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Table 3. Atomic data and solar abundances for the investigated 
lines. 



Ion 



A, A 



log gf Ref logCft Ref logEo 



Nai 


5889.96 


0.00 


0.11 


2 


-31.60 


14 


6.28 




5895.93 


0.00 


-0.19 


2 


-31.60 


14 


6.28 


Mgi 


4571.10' 


0.00 


-5.62 


2 


-31.96 


14 


7.58 




4702.99' 


4.33 


-0.44 


2 


-29.71 


14 


7.49 




5528.41' 


4.33 


-0.50 


2 


-30.20 


14 


7.53 




5711.09' 


4.33 


-1.72 


2 


-29.89 


14 


7.55 




5172.68 


2.70 


-0.45 


3 


-30.88 


14 


7.56 




5183.60 


2.70 


-0.24 


3 


-30.88 


14 


7.57 


Ah 


3961.53 


0.01 


-0.34 


2 


-31.20 


14 


6.44 


Ki 


7664.92 


0.00 


0.13 


4 


-31.00 


5 


5.12 




7698.98 


0.00 


-0.17 


4 


-31.00 


5 


5.12 


Cat 


4425.44 


1.87 


-0.36 


6 


-30.90 


* 


6.41 




4578.56 


2.51 


-0.70 


7 


-30.30 


16 


6.37 




5261.71 


2.51 


-0.58 


7 


-30.86 


16 


6.44 




5349.47 


2.70 


-0.31 


7 


-31.45 


16 


6.41 




5588.76 


2.51 


0.36 


7 


-31.39 


16 


6.28 




5590.12 


2.51 


-0.57 


7 


-31.39 


16 


6.39 




5857.45 


2.92 


0.24 


7 


-30.61 


16 


6.36 




6122.22 


1.88 


-0.32 


6 


-30.30 


15 


6.26 




6162.17 


1.89 


-0.09 


6 


-30.30 


15 


6.28 




6166.44 


2.51 


-1.14 


7 


-30.48 


16 


6.42 




6439.07 


2.51 


0.39 


7 


-31.58 


16 


6.27 




6471.66 


2.51 


-0.69 


7 


-31.58 


16 


6.34 




6493.78 


2.51 


-0.11 


7 


-31.58 


16 


6.31 




6499.65 


2.51 


-0.82 


7 


-31.58 


16 


6.38 




6449.81 


2.51 


-0.50 


7 


-31.45 


16 


6.34 


Can 


3933.66 


0.00 


0.10 


8 


-31.72 


17 


6.35 




8498.02' 


1.69 


-1.42 


8 


-31.51 


17 


6.35 


Een 


4508.29 


2.84 


-2.32 


9 


-32.00 


18 


7.44 




4582.84 


2.83 


-3.14 


9 


-32.03 


18 


7.43 




4923.93 


2.88 


-1.36 


9 


-32.05 


18 


7.41 




5018.44 


2.88 


-1.23 


9 


-32.05 


18 


7.45 




5197.57 


3.22 


-2.24 


9 


-32.05 


18 


7.46 




5234.62 


3.21 


-2.14 


9 


-32.05 


18 


7.41 




5264.81 


3.22 


-3.13 


13 


-32.03 


18 


7.53 




5325.56 


3.21 


-3.26 


12 


-32.05 


18 


7.56 




5425.25 


3.19 


-3.30 


12 


-32.05 


18 


7.51 




6247.56 


3.87 


-2.43 


13 


-32.00 


18 


7.57 




6456.38 


3.89 


-2.18 


13 


-32.00 


18 


7.59 


Srii 


4215.53' 


0.00 


-0.17 


10 


-31.80 


* 


2.92 


Ban 


4554.03' 


0.00 


0.16 


10 


-31.65 


* 


2.21 




6496.90 


0.60 


-0.38 


10 


-31.28 


17 


2.21 


Euii 


4129.72' 


0.00 


0.22 


11 


-32.08 


19 


0.56 



Excitation energy E„r is given in eV; 

' IS and HES components are taken into account; 

Ref.: 2 - NIST database; 3 - Aldenius et al. ([2007}; 4 - Butler ( [2000t ; 

5 - Zhang et al. ( I2006ab ; 6 - Smith & O'Neil ( fT975l ); 7 - Smith & 

Raggett ( fT98n i; 8 - Theodosiou ( fT989] l; 9 - Landstreet (fT998ll; 10 - 

Reader et al. fT980t ; 1 1 - Lawler et al. ( 120011 1; 12 - Moity ([1983 1; 13 

- Raassen & Uylings (1998); 14 - Gehren et al. ( [20041 ; 15 - Anstee & 

O'Mara {19951 1; 16 - Smith cl98U ; 17 - Barklem & O'Mara Cl998j ; 

18 - Barklem & Aspelund-Johansson {2005l >; 19 - Kurucz ( [T992] >; * - 

solar line profile fitting. 
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Table 4. Stellar parameters and element abundances of the se- 
lected stars where V„^ic is given in kms"' 



Table 6. Effect on element abundances (in dex) of BD -(-3°740 
and BD -13°3442 caused by uncertainties of its stellar param- 
eters 



HD/BD 


122563 
4600 


84937 
6365 


+3°740 
6340 


-13°3442 
6390 










T,ff,K 


Chemical 


reff(K) 


logs 


V„,„(kms-') 


log^ 


1.50±o.2 

-2.53±o.02 


4.00±0.12 

-2.15±0.04 


3.90±0.15 
-2.65±o.oi 


3.88±o.i5 
-2.66±o.oi 


species 


-100 


+0.16 -0.28 


+0.2 


[Fe/H] 


Fen 


-0.02 


+0.05 -0.09 


<0.01 


V ■ 

^ mic 


1.9 
NLTE 


1.6 


1.4 


1.4 


Cai 
Can 


-0.06 


+0.02 +0.01 


<0.01 


mm\ 


-2.22±0.06 


-1.76±0.04 


-2.35±o.03 


-2.26±0.05 


(^3933) 


-0.10 


-0.03 +0.04 


<0.01 


[Na/H] 


-2.93±0.0l 


-2.48±0.02 


-3.22±o.03 


-2.82±0.01 


(^8498) 


-0.08 


+0.06 -0.05 


-0.06 


[Al/H] 


-3.09 


-2.21 




-2.61 










[Cal/H] 


-2.32±0.06 


-1.72±0.04 


-2.19±0.05 


-2.12±0.06 










[Call/H] 


















(^3933) 


-2.42 


-1.70 


-2.16 


-2.15 










(.18498) 


-2.36 


-1.76 


-2.14 


-2.06 










[K/H] 


-2.60±o.04 


-1.87±0.05 


-2.37±0.04 


-2.28 










[Ba/H] 


-3.54 


-2.15 


2.85 


-2.89 










[Sr/H] 


-2.93 


-2.18 


-2.62 


-2.35 










[Eu/H] 


-3.04 
LTE 


-1.45 














[Mg/H] 


-2.39±o.io 


-1.88±o.03 


-2.40±0.05 


-2.29±0.03 










[Na/H] 


-2.51 ±0.04 


-2.02±0.03 


-3.03±0.02 


-2.58±o.08 










[Al/H] 


-3.13 


-3.04 




-3.47 










[Cal/H] 


-2.52±o.io 


-1.76±0.04 


-2.32±0.07 


-2.27±0.06 










[Call/H] 


















(^3933) 


-2.42 


-1.70 


-2.16 


-2.15 










(/18498) 


-2.29 


-1.56 


-1.92 


-1.81 










[K/H] 


-2.25±0.02 


-1.65±o.03 


-2.17±0.04 


-2.08 










[Ba/H] 


-3.57 


-2.30 


-3.32 


-3.39 










[Sr/H] 


-2.85 


-2.18 


-2.89 


-2.57 










[Eu/H] 


-3.16 


-1.61 















Table 5. Effect on element abundances (in dex) and the Ba odd 
isotope fraction fodd (in %) of HD 122563 caused by uncertain- 
ties of its stellar parameters. 



Chemical 


reff(K) 


lo 


%S 


V„,„(kms-') 


species 


+ 120 


-0.2 


+0.2 


+0.1 


Fen 


< +0.01 


-0.07 


+0.08 


<0.01 


Cai 


+0.11 


+0.03 


-0.03 


<0.01 


Can 










(^3933) 


+0.11 


+0.02 


-0.01 


<0.01 


(^8498) 


+0.12 


+0.00 


+0.09 


<0.01 


Mgi 


+0.11 


+0.03 


-0.03 


-0.02 


Ban 


+0.07 


-0.08 


+0.08 


<0.01 


fodd 


+ 12 


+2 


-3 


-8 



